Non-equilibrium quantum condensation in an incoherently pumped dissipative system 
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We study spontaneous quantum coherence in an out of equilibrium system, coupled to multiple 
baths describing pumping and decay. For a range of parameters describing coupling to, and occu- 
pation of the baths, a stable steady-state condensed solution exists. The presence of pumping and 
decay significantly modifies the spectra of phase fluctuations, leading to correlation functions that 
differ both from an isolated condensate and from a laser. 
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The phenomenon of condensation, i.e macroscopic oc- 
cupation of a single quantum mode, has attracted con- 
siderable attention in recent years. It ranges from Bosc- 
Einstcin Condensation (BEC) of structureless bosons to 
the BCS-type collective state of fermions and has been 
studied in several physical systems such as degenerate 
atomic gases and superconductors pj. Further, recent 
experimental advances in manipulation of atomic Fermi 
gases have led to realisation of the BCS-BEC crossover 
regime 0. The next challenge is to control and study 
condensed states in solid-state. For this, the currently 
promising candidates arc 1 excitons in coupled quantum 
wells , microcavity polaritons 0, 0, E| , quantum Hall 
bilayers hj, and Josephson junction arrays in microwave 
cavities yj- 

Unlike atomic gases, solid-state systems face dephas- 
ing and decay, as (with the special exception of supercon- 
ductors) it is not usually possible to isolate the conden- 
sate from the environment: P honoris and impurities lead 
to dephasing, and due to imperfect trapping, particles 
escape, requiring external pumping to sustain a steady- 
state. If such processes are faster than thermalisation the 
system remains out of thermal equilibrium. Dissipation 
and decay not only present experimental obstacles, but 
also pose fundamental questions about the robustness of 
a condensate: Is a steady-state condensate possible with 
incoherent pumping and decay, if so, how does it differ 
from thermal equilibrium? Condensation in dissipative 
systems also provides a connection to the laser 0. The 
relation between lasing and BEC is particularly relevant 
for polariton BEC, where the experimental distinction 
between the two is not straightforward [lfij ]. 

Models that combine potentially strong non- 
equilibrium pumping with spontaneous symmetry 
breaking are not well-studied, which motivates a study 
of simple models to extract the principal qualitative 
features. In this Letter we study spontaneous con- 
densation in a system coupled to independent baths, 
not in thermal or chemical equilibria with each other, 
providing incoherent pumping and decay. We focus on a 
Bose-Fermi system with disorder localised fermions; this 
is a model for exciton-polaritons 0, 0, or Josephson 



junctions in microcavities However, many of our 
conclusions apply more generally to condensation with 
pumping and decay. We show that steady-state sponta- 
neous condensation can occur in such systems, and can 
be distinct from lasing: The condensate can exist at low 
densities, far from the inversion required for lasing. We 
study fluctuations about a steady-state condensate and 
find that the collective modes are qualitatively altered 
by the presence of pumping and decay: The low energy 
phase mode (Goldstone, Bogoliubov mode) becomes 
diffusive at small momenta. By considering the effect 
of phase fluctuations, we find the decay of correlations, 
which at large times and distances differs both from that 
for a thermal equilibrium condensate and from a laser. 
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describes two fermionic species b a and a ai interacting 
with bosonic modes ip p , and coupled to three baths; 
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Efc ^k^k^k- A single two level system coupled to multi- 
ple baths has been recently studied in the context of the 
Kondo problem . Condensed solutions of eq. ([Til have 
been studied in the context of atomic Fermi gases [12| and 
microcavity polaritons 0, Q| . In this Letter we focus 
on microcavity polaritons, so bl t ,a a describe an electron 
and hole within a disorder-localised exciton state of en- 
ergy e a . This can also be viewed as a fermionic represen- 
tation of a hard-core boson or of a spin. These arc dipole 
coupled to cavity photon modes ip p , with low p disper- 
sion, lo p ~ loq + p 2 /2m p h, where m p h = (h/c)(2w/w) is 
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the photon mass in a 2D microcavity of width w. Duo 
to the finite reflectivity of the cavity mirrors, photons es- 
cape, so the system must be pumped (excitons injected) 
to sustain a steady-state. Incoherent fermionic pump- 
ing and photon decay are described by where A k , B k 
arc fermionic annihilation operators for the pump baths, 
while $ k are bosonic annihilation operators for photon 
modes outside the cavity. 

With pumping and decay i.e. coupling to baths with 
different temperatures and/or chemical potentials, the 
system distribution function can be far from thermal and 
needs to be obtained self-consistently with the system's 
spectrum. Therefore we use non-equilibrium field theory 
in the Keldysh path-integral formulation Integrat- 
ing over the bath's degrees of freedom and fermion fields 
yields an effective action for photons. The dependence 
of the effective action on the bath properties is parame- 
terised by the functions k(u>), Fa,b(u), and Fy(uj). 

The cavity decay rate k(lu) = tv£ 2 (u))N''(<jj), where £ is 
the coupling of the cavity photons to the bosonic modes 
in eq. (0 and N 1 * is the density of states of these modes. 
Similarly 7(0;) = nT 2 (uj)N r (ui) where T and N r arc re- 
spectively coupling to, and the density of states of, the 
fermionic pumping baths. In this work we assume a flat 
spectrum for the baths, so k and 7 are frequency inde- 
pendent. Frequency dependence is however present in the 
bath distribution functions; Fa,b(&) = 1 — 2n A ' B (ui), and 
Fij,(u)) = I + 2n*(o>), where n A ' B ,71* are occupations of 
the baths. 

We proceed by extremising the non-equilibrium effec- 
tive action with a steady-state, uniform, photon field of 
the form, ip(t) = ipe~ %tJ,st 1 giving a non-equilibrium gen- 
eralisation of the usual gap equation (e.g. fujlH^: 



(u - fig - iK)i/j = gTr(iGg a ) 
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Here, iGjf a (t,t') = (a cl (t)b cl (t')), is the anomalous 
Keldysh fermionic Green's function, where a c i = (a/ + 
a&)/ \/2 and /, b are the forward and backward branches 
of the Keldysh time contour and is given by: 
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(F A + F b )uj + {F b - F A )(e + ij) 
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where E = \Je 2 + g 2 \ip\ 2 : e = e — fig/2, and the ar- 
guments of Fb,a(lo) are shifted by ±)ig/2 while lo is a 
real frequency. The bosonic bath's distribution Fy does 
not enter the mean-field gap equation, as the mean-field 
considers only the condensed photons, and not the pop- 
ulation of incoherent photons. In the limit 7, k — > with 
the bath distributions Fa,b being thermal, eq. (jSJ) re- 
duces to its equilibrium form 0,13b f° r finite 7 and k 
it is significantly altered. 

As in thermal equilibrium, the normal state ip = is 
always a solution of eq. ©, but for some range of pa- 
rameters there is also a condensed ip 7^ solution. When 
a solution ip ^ exists, the solution ip ~ becomes 



unstable. To understand this instability, we consider 
small fluctuations about the mean-field, ip — ipo + Sip. 
The effective action for these fluctuations Sip has a part 
from the free photon action, and a part from interactions 
with fermions, ■^Tr(G5G~ 1 G5G~ 1 ), where the fermionic 
Green's functions G are four by four matrices in the 
Keldysh and particle-hole (a, b) spaces. Inverting the 
effective action for fluctuations gives the photon Green's 
functions, J) K ' R ' A . The retarded and advanced Green's 
functions V R ' A give the excitation spectrum. The distri- 
bution function F s , defined by V K = V R F S - F S V A , 
determines how the spectrum is occupied. As the system 
need not be in thermal or chemical equilibria with the 
baths, Fg is in general not thermal, and differs from the 
bath distributions, Fa,b and Fq,. 

The instability of the normal state when a condensed 
solutions exists is analogous to that in thermal equi- 
librium. Even when the bosonic distribution Fg is 
far from thermal, as the system approaches the phase 
transition, Fg will have a divergence at a frequency 
which we define as an effective chemical potential, where 

SOD^Vcff)) = °- Taking the zeros of 3?(P fl_1 (w*, q)) 
as defining the normal modes of the system, condensa- 
tion occurs when /i e g reaches the bottom of the band 
of excitations, as in equilibrium. This condition, that a 
solution to D R (M e ff , q = 0) = exists, is equivalent to 
the gap equation, eq. ©, at tp = 0, /i c g- = fig, since D R 
describes a susceptibility which diverges at the transi- 
tion. Beyond this point, the normal state is unstable, as 
/i e g? would lie in a bosonic band. This idea of instability 
can be directly connected to another definition: Beyond 
this point, the poles of the D R have positive imaginary 
parts, fluctuations grow (rather than decay) in time ex- 
ponentially. To see this, consider the imaginary parts of 
the poles of T> as a function of momentum, q. At large 
q, these poles describe bare photons, and so are stable. 
By the previous definitions, above the transition, there is 
a q at which the normal state Green's function has a real 
pole, V R (lu* , q) = 0. The sign of the imaginary part of 
the pole changes at this point, so the low q poles are un- 
stable. Thus, whenever there is a condensed solution to 
the gap equation the ip — solution is always unstable. 

There is, however, an important difference between the 
non-equilibrium steady-state and thermal equilibrium. 
Unlike thermal equilibrium, there is a range of param- 
eters for which neither the normal state nor condensed 
solutions of the form ip(t) — ipe~ l ^ st : are stable. This 
is not too surprising as systems similar to JQ) are known 
to follow a complicated or even chaotic dynamics [l7j |. 
Although such anharmonic solutions would be of a great 
interest, we focus here on steady-state condensed solu- 
tions of the usual form. 

Having understood the stability of solutions, we next 
solve the gap equation Q to find the non-equilibrium 
phase diagram. Since eq. © is complex it gives two 
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equations for two unknowns: the order parameter tjj and 
the frequency ps- The common oscillation frequency ps 
would in thermal equilibrium be the system's chemical 
potential, considered as a control parameter, adjusted to 
match the required density, and the (real) gap equation 
determines only -0. Here, because different baths have 
different chemical potentials, the system is not in chemi- 
cal equilibrium with either bath, so both ps and ip must 
be found from the gap equation. The density, which can 
be found given ip and ps, is set by the relative strength 
of the pump and decay 

From here, the bath distributions are taken to be indi- 
vidually in thermal equilibrium, although the formalism 
would allow any distribution. However, as the baths need 
not be in equilibrium with each other, the system can still 
be far from thermal equilibrium. Further, for simplicity, 
the figures shown are all for the baths at zero tempera- 
ture, so the bath distributions are defined by their chem- 
ical potentials. To ensure that on average only one of the 
two fermionic levels a, b is occupied, the chemical poten- 
tials of baths A, B are related by pa = —Pb- Since the 
cavity photon modes start at energies much above the 
zero for bulk photon modes, we take the chemical poten- 
tial of the decay bath to be large and negative. With 
these restrictions, the baths are described by three pa- 
rameters, the couplings 7, k and ps parameterising the 
occupation of the pumping baths. 

Figure^shows a phase diagram in terms of the param- 
eters 7, k, fiB- It shows both where a condensed solution 
to the gap equation exists, and where the solution is sta- 
ble. At a given decay rate k there is a minimum 7 (as 
pumping is proportional to 7) and a maximum 7 (as de- 
phasing is also proportional to 7) required for condensa- 
tion. Stable condensed solutions exist only for k smaller 
than about Q.2g. In the region shown in Fig. ^ the con- 
densed solutions are all below population inversion. How- 
ever for 7 > g, when system is in a weak coupling regime, 
only laser-like solutions which require population inver- 
sion are possible. For large 7 and pb (large pumping), 
our theory recovers the regular laser limit. 

We now discuss the collective modes, considering fluc- 
tuations about the steady-state. Motivated by mi- 
crocavity polaritons, we study these collective modes 
by calculating the photoluminescence (PL) spectra: 
W^,(t,t') — (fpf(t)tp b (t')). In the normal state, as ex- 
pected, one finds homogeneously broadened upper and 
lower polariton modes. This broadening depends on 
all of the bath parameters, k,7 and ps- Approaching 
the phase boundary from the normal side, the lower 
polariton lincwidth reduces to zero. This can be un- 
derstood by identifying the zeros of the real part of 
the inverse Green's function as the polariton energies, 

the lincwidth, 



0, and the imaginary part as giving 
1/t p w 3(X> fl_1 (w*,g)). The earlier dis- 
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FIG. 1: Mean field phase diagrams. Left: Critical 7 vs decay 
rate «, for two different pumping bath chemical potentials, 
fi b ■ Right: Critical 7 vs ps for K=0.038g as in experiment [jj. 
Solid lines show where a solution of the gap equation exists, 
dashed lines show where that solution is stable. The region 
where a stable condensed solution exists is shaded. Stars mark 
the choice of parameters shown in Fig. [5] 



implies that 1/t p vanishes at the transition. 

In the condensed state, as in equilibrium^!], the ex- 
citation spectrum changes. This leads to a soft mode, 
describing phase fluctuations, as global phase rotation 
symmetry is broken. These phase fluctuations may be 
large, so one must include the phase fluctuations to all 
orders 0] to find the field-field correlations; amplitude 
fluctuations remain small, as they have a restoring force. 
Writing, ip(t) = ^ p Q + 7r(i)e 11 ^*), where p is the mean- 
field condensate density, the PL spectrum is thus: 

OW*> r) = P0 [l + O (1/po)] cxp [-f(t, r)] , (5) 



f(t,r)= / dv (dqf 1 



AiV<{v,q). (6) 



The phase-phase Green's function X>^ is found by in- 
verting the action expanded to quadratic order. 

The 0(1/ po) terms in eq. JjjJ, due to amplitude- 
amplitude and phase-amplitude Green's functions, are 
included when plotting Fig. [3 The form of these terms, 
and the full expression for 27^,, are complex, so we do 
not reproduce them here. The behaviour of 2?^ in the 
limit of small frequencies and momenta, is of interest, 
and takes a simple form: 



Av 2 x 2 



(7) 



cussion relating the gap equation to T> R (ps,<l = 0) form v = —ix±i^Jx 



where C, c, and x can be found from the full expressions. 
Without pumping and decay, the term x would vanish. 
Its presence means that rather than a linear dispersion 
at low q, as in a closed system in thermal equilibrium^^ 
llfij . the poles of the Green's function have the following 

For \q\ <x/c, these modes 



c 2 q 2 . 
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are diffusive; the poles are entirely imaginary, only for 
\q\ > x/c, do they acquire a real part. This behaviour 
is apparent in the PL spectrum in Fig. [21 The diffusive 
behaviour at small q is not limited to Bose-Fermi systems; 
it should also be present in other condensed systems with 
pumping and decay. Note that at q = and v = 0, 2?^ 
has a real pole; a manifestation of broken symmetry in 
the infinite system. 

As given by eq. (jHJ), the PL spectrum does not de- 
pend linearly on T>^. If 2?^ were small, the exponen- 
tial in eq. (J5J could be expanded, giving PL divided be- 
tween a condensate term po5(y)5(q), and a part from 
occupation of the fluctuation spectrum. However, at 
v, q — ► 0, 2?00 is not small, and so phase fluctuations 
give a lineshape to the condensate and determine long 
time field-field correlations. In 2D, inserting eq. Q in 
eq. ©, the main dependence of f(t,r) on t,r comes 
from a logarithmic divergence, J dq/q, cut at large mo- 
menta by terms beyond those in eq. JjJ, and at small 
momenta by one of 1/r, 1/ct , 1 jcJtT g or \/cr g , where 
T g = 1/x is the lifetime of the phase mode. This loga- 
rithmic form leads to power law field-field correlations. 
(This is true also in thermal-equilibrium^^, for which, 
f{t,r) = rjh-i(Vc 2 t 2 + r 2 //3c) , with rj oc k B T/p .) Ac- 
cording to the relative values of r,t,T g , different lower 
cutoffs apply, and so the power law differs at different 
places in the r, t plane. For small r, t we recover the equi- 
librium power laws; however when f>T 9 , and c^JtT g 3> 
r, then r g becomes important, and f(t) oc \a{cy/Urg~) , giv- 
ing a condensate lineshape that differs from that of an 
non-dissipative, thermal equilibrium 2D system. Power 
law field correlations at long times lead to a power law 
divergence of their Fourier transform as u, q — » 0, so there 
is no well defined linewidth. This differs from phase dif- 
fusion for a single mode studied in Laser theory |<J], where 
phase correlations grow linearly in time, giving exponen- 
tial decay of field correlations. 

At higher q and v (FigEllnsets) the difference between 
expanding to linear order in Green's functions and keep- 
ing all orders of phase fluctuations is not visible. The 
strong coupling spectrum consists of the usual phase and 
amplitude modes (marked a,b) emerging from the lower 
and upper polariton branches. One can also see an oc- 
cupation edge (c), and above that the lower polariton 
following the exciton dispersion (d). These features at 
higher momenta are analogous to thermal equilibrium 
non-dissipative polariton condensation in strong coupling 
regime ll-il but with non-thermal occupations. Increas- 
ing pumping, dephasing increases, and the system crosses 
to weak-coupling (lower panels of Fig. El, the polariton 
splitting is suppressed, and the spectrum follows the pho- 
ton dispersion. The diffusive region, then linear disper- 
sion, at small momentum, seen in the main figures, oc- 
curs both in strong and weak coupling and is a sign of 
condensation in the dissipative system. 

To conclude, we have studied how steady-state sponta- 
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FIG. 2: Photoluminescence iT>^ Av,q) of a condensed sys- 
tem, where q is shown by angle of emission tan -1 (cq/wo). 
Top: Strong coupling, Bottom: weak coupling, exact param- 
eters marked by stars on Fig. Main figures show PL from 
small v, q region, to all orders in phase fluctuations. Insets 
show a larger range of v, q for the same parameters. (Dotted 
lines have been added to show the faint amplitude mode). 

neous condensation emerges in non-equilibrium systems 
with pumping and decay. This condensation is distinct 
from lasing: It can occur at densities much lower than the 
population inversion required for lasing, and the decay of 
correlations, and thus condensate lineshape, differ from 
that for phase diffusion of a single laser mode. Dissipa- 
tion qualitatively changes the fluctuation spectrum with 
respect to isolated thermal equilibrium: At low momenta, 
the phase mode becomes diffusive, changing the power- 
laws controlling long time decay of field-field correlations. 
These conclusions, although studied here for microcavity 
polaritons, apply also to other Bose and Bose-Fermi con- 
densates subject to pumping and decay. 
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